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Abstract. We investigate the phase diagram of the so-called Polyakov-Nambu-Jona-Lasinio model at finite temperature 
and nonzero chemical potential with three quark flavours. Chiral and deconfinement phase transitions are discussed, and the 
relevant order-like parameters are analyzed. The results are compared with simple thermodynamic expectations and lattice 
data. A special attention is payed to the critical end point: as the strength of the flavour-mixing interaction becomes weaker, 
the critical end point moves to low temperatures and can even disappear. 
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There is strong evidence that quantum chromodynam- 
ics (QCD) is the fundamental theory of strong interac- 
tions. Its basic constituents are quarks and gluons that 
are confined in hadronic matter. It is believed that at high 
temperatures and densities hadronic matter should un- 
dergo a phase transition into a new state of matter, the 
quark-gluon plasma (QGP). A challenge of theoretical 
studies based on QCD is to predict the equation of state, 
the critical end point (CEP) and the nature of the phase 
transition. 

As the equation of state of strong interacting matter 
given by lattice simulations is now considered as a func- 
tion of temperature and a limited range of chemical po- 
tential, QCD-like models, as Nambu-Jona-Lasinio (NJL) 
type models, are suited to provide guidance and informa- 
tion relevant to observe experimental signs of deconfine- 
ment and QGP features. 

NJL type models, that take into account only quark de- 
grees of freedom, give the correct chiral properties; static 
gluonic degrees of freedom are then introduced in the 
NJL Lagrangian through an effective gluon potential in 
terms of Poly akov loop with the aim of including features 
of both chiral symmetry breaking and deconfinement. 

Hence our calculations are performed in the frame- 
work of an extended SU(3)/ PNJL Lagrangian, which 
includes the 't Hooft instanton induced interaction term 
that breaks the Ua(1) symmetry, and the quarks are 
coupled to the (spatially constant) temporal background 
gauge field <J> [1, 2]. The Lagrangian reads: 
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The covariant derivative is defined as D 11 = d 11 — iA^, 
with A' 1 = 8qAq (Polyakov gauge); in Euclidean notation 
Aq = —1A4. The strong coupling constant g is absorbed 
in the definition of A'' (x) = gsrfa (x)^f, where s/a is the 
(SU(3) C ) gauge field and A„ are the (color) Gell-Mann 
matrices. 

The effective potential for the (complex) field <t> 
adopted in our parametrization of the PNJL model reads: 
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The effective potential exhibits the feature of a phase 
transition from color confinement (T < Tq, the minimum 
of the effective potential being at <t> = 0) to color decon- 
finement (T > Tq, the minima of the effective potential 
occurring at <t> ^ 0). 

The parameters of the effective potential are given 
by oo = 3.51, a\ = —2.47, 02 = 15.2 and b^ = —1.75. 
These parameters have been fixed in order to reproduce 
the lattice data for the expectation value of the Polyakov 
loop and QCD thermodynamics in the pure gauge sector. 
The parameter Tq, the critical temperature for the decon- 
finement phase transition within a pure gauge approach, 
was fixed to 270 MeV, according to lattice findings. This 
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FIGURE 1. Scaled pressure (p), energy per particle (e), and entropy (s) as a function of the temperature at zero 
chemical potential. The data points are taken from [4]. 



choice ensures an almost exact coincidence between chi- 
ral crossover and deconfinement at zero chemical poten- 
tial, as observed in lattice calculations. We notice, how- 
ever, that a rescaling of Ty to 210 MeV may be needed in 
some cases in order to get agreement between model cal- 
culations and thermodynamic quantities obtained on the 
lattice (see Fig. 1). Let us stress that this modification of 
Tq is essentially done for rescaling reasons (an absolute 
temperature scale has not a very strong meaning in these 
kind of Ginzburg-Landau model based on symmetry) but 
does not change drastically the physics. 

The parameters of the NJL sector are: m u = m £ / = 
5.5 MeV, m s = 140.7 MeV, 8s A 2 = 3.67, gD A 5 = - 12.36 
and A = 602.3 MeV, which are fixed to reproduce the 
values of the coupling constant of the pion, f n = 92.4 
MeV, and the meson masses of the pion, the kaon, the rj 
and tj', respectively, M K = 135 MeV, M K = 497.7 MeV, 
Mn = 514.8 MeV andM,,, = 960.8 MeV. 

In the limit of vanishing quark chemical potential, 
significant information on the phase structure of QCD 
at high temperature is obtained from lattice calculations. 

In Fig. 1, we plot the scaled pressure, the energy and 
the entropy as functions of the temperature compared 
with recent lattice results (see Ref. [4]). Since the tran- 
sition to the high temperature phase is a rapid crossover 
rather than a phase transition, the pressure, the entropy 
and the energy densities are continuous functions of the 
temperature. We observe a similar behavior in the three 
curves: a sharp increase in the vicinity of the transition 
temperature. 

Asymptotically, the QCD pressure for Nf massless 
quarks and (N 2 — 1 ) massless gluons is given (p. = 0) 
by: 

PM =(N 2_ l) El +NcN ^ (4) 
T 4 y c ; 45 ; 180 

where the first term denotes the gluonic contribution and 
the second term the fermionic one. The results follow 
the expected tendency and go to the free gas values 
(or Stefan-Boltzmann (SB) limit). The inclusion of the 



Polyakov loop effective potential ( < t , ,<i>; T) (it can be 
seen as an effective pressure term mimicking the gluonic 
degrees of freedom of QCD) is required to get the correct 
limit (indeed in the NJL model the ideal gas limit is 
far to be reached due to the lack of gluonic degrees of 
freedom). 

Some comments are in order concerning the role of 
our regularization procedure for T > T c , that allows for 
the presence of high momentum quark states [3]. In this 
temperature range, due to the presence of such states, 
the physical situation is dominated by the significant 
decrease of the constituent quark masses. This allows for 
an ideal gas behavior of almost massless quarks with the 
correct number of degrees of freedom. 

Let us notice that, just below T c , the pressure and the 
energy fail to reproduce the lattice points: for example 
there is a small underestimation of the pressure and 
energy in the model calculations. It was pointed out that 
the lack of mesonic correlations in the PNJL model is 
responsible for, at least, a fraction of this discrepancy. 

Our numerical results for the isentropic lines in the 
(T,n) plane are shown in Fig. 2 (left side). We start the 
discussion by analyzing the behavior of the isentropic 
lines in the limit T — > 0. Our convenient choice of the 
model parameters allows a better description of the first 
order transition than other treatments of the NJL (PNJL) 
model. This choice is crucial to obtain important results: 
the criterion of stability of the quark droplets is fulfilled, 
and, in addition, simple thermodynamic expectations in 
the limit T — > are verified. In fact, in this limit s — > 
according to the third law of thermodynamics and, as 
Pb —> too, the satisfaction of the condition s/ps = 
const, is insured. 

We observe that the isentropic lines with s/pb = 
1,...,6 come from the region of symmetry partially re- 
stored and attain directly the phase transition. Conse- 
quently, all isentropic trajectories directly terminate in 
the end of the first order transition line at T = 0. The 
trajectories with s/ps > 6 go directly to the crossover re- 




FIGURE 2. Left: Isentropic trajectories in the (T, fl) plane. Right: Dependence of the location of the CEP on the strength of the 
t Hooft coupling constant go- 



gion and display a smooth behavior, although those that 
pass in the neighborhood of the CEP show a slightly kink 
behavior [3]. 

The location and even the existence of the CEP in the 
phase diagram is a matter of debate. While different lat- 
tice calculations predict the existence of a CEP [5], the 
absence of the CEP in the phase diagram was seen in re- 
cent lattice QCD results [6], where the first order phase 
transition region near ji = shrinks in the quark mass 
and ji space when ji is increased [6]. Due to the im- 
portance of the Ua(1) anomaly and its influence on sev- 
eral observables, it is demanding to investigate possible 
changes in the location of the CEP in the (T, ju) plane 
when the anomaly strength is modified. In Fig. 2 (right 
side) we show the location of the CEP for several values 
of gD compared to the results for go , the value used for 
the vacuum. As already pointed out by K. Fukushima in 
[1], we also observe that the location of the CEP depends 
on the value of go- In fact, our results show that the ex- 
istence or not of the CEP is determined by the strength 
of the anomaly coupling, the CEP getting closer to the jj. 
axis as go decreases. 

We investigated the phase diagram of the so-called 
PNJL model at finite temperature and nonzero chemi- 
cal potential with three quark flavours. Chiral and de- 
confinement phase transitions are discussed, and the rel- 
evant order-like parameters are analyzed. The results are 
compared with simple thermodynamic expectations and 
lattice data. A special attention is payed to the critical 
end point: as the strength of the flavour-mixing interac- 
tion becomes weaker, the critical end point moves to low 
temperatures and can even disappear. 

The sets of parameters used is compatible with the for- 
mation of stable droplets at zero temperature, insuring 
the satisfaction of important thermodynamic expectation 



like the Nernst principle. Consequently, all the trajecto- 
ries directly terminate in the same point of the horizontal 
axes at T — 0. The picture provided here is a natural re- 
sult in these type of quark models with no change in the 
number of degrees of freedom of the system in the two 
phases. 

Other important role is played by the regularization 
procedure which, by allowing high momentum quark 
states, is essential to obtain the required increase of ex- 
tensive thermodynamic quantities, insuring the conver- 
gence to the SB limit of QCD. In this context the gluonic 
degrees of freedom also play a special role. 

The successful comparison with lattice results shows 
that the model calculation provides a convenient tool 
to obtain information for systems from zero to nonzero 
chemical potential which is of particular importance to 
the knowledge of the equation of state of hot and dense 
matter. 
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